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Anthrax  protective  antigen  (PA)  is  one  of  the  three  proteins  produced  by  the  gram  positive 
bacteria  Bacillus  anthracis  collectively  known  as  the  “anthrax  toxin”  (Ascenzi,  P.;  Visca,  P.; 
Ippolito,  G.;  Spallarossa,  A.;  Bolognesi,  M.;  et  al.  Anthrax  toxin:  a  tripartite  lethal  combination. 
FEBS  Lett .  2002,  531 ,  384-388).  The  role  played  by  PA  in  anthrax  intoxication  is  to  transport 
the  two  enzymes  lethal  factor  (LF)  and  edema  factor  (EF)  into  the  cell.  Collier  and  co-workers 
(Mourez,  M.;  Kane,  R.  S.;  Mogridge,  J.;  Metallo,  S.;  Deschatelets,  P.;  et  al.  Designing  a  polyvalent 
inhibitor  of  anthrax  toxin.  Nat.  Biotechnol.  2001,  958).  reported  the  isolation  of  two  peptides 
via  phage  display  that  bind  to  the  PA63  heptamer  and  inhibit  its  interaction  with  LF  and  EF, 
and  thereby  prevent  the  transport  of  LF  and  EF  into  the  cell.  One  of  these  peptides,  His-Thr- 
Ser-Thr-Try-Trp-Trp-Leu-Asp-Gly-Ala-Pro  (PI),  was  selected  for  structural  investigation  on 
the  basis  of  its  ability  to  prevent  the  binding  of  LF  to  the  PA63  heptamer  bundle.  Two- 
dimensional  trNOESY  experiments  coupled  with  NOE  restrained  simulated  annealing  calcula¬ 
tions  were  used  to  determine  the  PA^-bound  conformation  of  PI.  On  binding  to  PAe3,  PI  adopts 
a  helical  conformation  involving  residues  3-9  while  the  C-  and  N-terminal  residues  exhibit 
dynamic  fraying. 


Introduction 

Anthrax,  from  the  Greek  word  for  coal  ( anthrakitis ), 
was  the  name  originally  given  to  the  cutaneous  illness 
caused  by  the  Gram-positive  bacteria  Bacillus  anthracis 
due  to  the  formation  of  black  skin  lesions.1  The  name 
now  applies  to  all  three  forms  of  anthrax:  cutaneous, 
gastrointestinal,  and  inhalational.  Most  toxins  that  act 
intracellularly,  such  as  ricin  and  the  botulinum  neuro¬ 
toxins  (BoNTs),  are  characterized  as  AB  toxin,  meaning 
that  the  toxin  molecule  contains  two  subunits,  or 
protomers,  that  are  released  after  proteolytic  cleavage 
of  the  holotoxin.1-3  Protomer  A  is  the  activating  or 
enzymatic  region,  and  B  is  the  binding  protomer,  which 
interacts  with  a  toxin-specific  receptor  located  on  the 
surface  of  the  target  cell.4  5  B.  anthracis  is  also  charac¬ 
terized  as  an  AB  toxin;  however,  instead  of  one  protein, 
the  bacterium  secretes  three  separate  proteins:  protec¬ 
tive  antigen  (PA:  83  kDa),  edema  factor  (EF:  89  kDa), 
and  lethal  factor  (LF:  83  kDa).  The  combined  actions 
of  these  proteins  constitute  the  anthrax  toxins  which 
induce  cell  death.5*6  It  should  be  pointed  out  that  each 
protein  by  itself  is  nontoxic.1  The  PA  serves  the  role  as 
the  common  binding  protomer  (protomer  B),  while  EF 
and  LF  serve  as  independent  A  protomers,  carrying  out 
different  enzymatic  reactions  within  the  host  cell.1*3 
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An  individual  PA  molecule  consists  of  four  structural 
domains  existing  mainly  as  antiparallel  ^-sheets.1*6 
Domain  1  contains  binding  sites  for  two  calcium  ions 
as  well  as  the  cleavage  site  for  the  furin  protease.1 
Domain  2  is  involved  in  pore  formation  and  contains  a 
large  flexible  loop  that  is  believed  to  be  involved  in 
membrane  insertion.1  Domain  3  is  involved  in  the 
formation  of  the  heptamer  bundle.  Last,  domain  4  binds 
to  the  anthrax  toxin  receptor  (ATR)  which  is  located  on 
the  surface  of  the  host  cell.1 

The  first  step  in  anthrax  intoxication  is  the  binding 
of  PA  to  the  ATR.  At  the  receptor,  PA  undergoes  a 
proteolytic  cleavage  by  furin,  or  a  furin-like  protease, 
at  a  surface  loop  within  domain  1  to  yield  two  frag¬ 
ments,6  a  20  kDa  N-terminal  fragment  (PA20)  and  a  63 
kDa  C-terminal  fragment  (PA^).  PA$3  remains  bound 
to  the  receptor  and  undergoes  self-assembly  with  other 
PA63  molecules  forming  a  ring-shaped  heptamer  bundle 
that  binds  one  to  three  molecules  of  EF  or  LF.  The 
interaction  of  the  seven  PA63  molecules  has  been 
described  as  packing  “pie  wedges”1  with  the  remainder 
of  domain  1  (after  proteolytic  cleavage)  and  domain  2 
lining  the  inside  of  the  bundle.  Domains  3  and  4  are 
aligned  on  the  outside  of  the  heptamer  bundle.1  This 
arrangement  presents  a  large  hydrophobic  region  for  the 
binding  of  LF  and  EF.1 

The  resulting  complex  is  transported  into  the  cell  via 
a  receptor-mediated  endocytosis,  and  is  subsequently 
carried  to  an  acidic  compartment  within  the  cell.1  The 
low  pH  causes  the  PA$3  heptamer/LF  (or  heptamer/EF) 
complex  to  cross  the  membrane  into  the  cytosol.  EF  is 
a  calcium-  and  CaM-dependent  adenylyl  cyclase  that 
results  in  an  increase  of  intracellular  cAMP,  which,  in 
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A)  B) 


C)  D) 


Figure  1.  The  !H  spectra  of  (A)  4.76  x  10  5  M  PA^,  (B)  4.76  x  10  5  M  PA6.3  plus  5.23  x  10  4  M  PI  (peptide  to  protein  ratio  11:1) 
(C)  4.76  x  10-5  M  PA63  plus  1.23  x  10~3  M  PI  (peptide  to  protein  ratio  25:1),  and  (D)  8.3  x  10“4  M  PI  (all  in  150  mM  sodium 
acetate  buffer,  pH  =  4.0,  277  K).  As  can  be  seen  in  this  figure  there  are  only  minor  detectable  changes  in  chemical  shift  on 
increasing  peptide  concentration.  A  comparison  of  spectrum  C  with  spectrum  D  shows  a  broadening  of  the  resonances.  This  is  a 
result  of  the  contribution  of  the  bound  state  T2  relaxation  rate  as  compared  to  the  T2  relaxation  of  the  free  state.14  41  The  observed 
chemical  shift  and  line  broadening  behavior  is  consistent  for  the  peptide  in  fast  exchange  between  a  bio-molecule-bound  and  free 
forms. 


turn,  causes  disruption  of  intracellular  communica¬ 
tion.7’8  EF  also  disrupts  neutrophil  function,  resulting 
in  edema.  In  general,  death  of  the  host  is  not  believed 
to  be  induced  by  the  action  of  EF.7’8  LF  is  a  metallo- 
protease  that  cleaves  mitogen-activated  protein  (MAP) 
kinase-kinases,  MEK1  and  MEK2,  at  their  N-terminus 
by  hydrolyzing  the  Pro8-Ile9  amide  bond  of  MEK1  and 
the  Pro10- Ala11  amide  bond  of  MEK2.910  This  most  likely 
results  in  the  overproduction  of  certain  lymphokines  by 
the  macrophages,  and  thereby  causes  cell  death  through 
a  type  of  septic  shock.6 

Collier  and  co-workers2  reported  the  isolation  of  two 
peptides  via  phage  display  that  bind  to  the  PA63 
heptamer  and  inhibit  its  interaction  with  LF  and  EF, 
and  thereby  prevent  the  transport  of  LF  and  EF  into 
the  cell.  The  two  peptides2  consist  of  the  following  amino 
acid  sequences:  His-Thr-Ser-Thr-Tyr-Trp-Trp-Leu-Asp- 
Gly- Ala-Pro  (PI)  and  His-Gln-Leu-Pro-Gln-Tyr-Tyr-Trp- 
Trp-Leu-Ser-Pro-Gly  (P2).  PI  was  reported  to  inhibit  the 
binding  of  radiolabeled  LFn  (a  255-amino  acid  residue 
peptide  corresponding  to  the  N-terminal  region  of  the 
PA{>3  binding  domain  of  LF)  to  PA63  on  Chinese  hamster 
ovary  (CHO)  cells  with  an  IC50  of  approximately  150 
/*M.2  Collier  and  Whitesides2  developed  inhibitors  of  the 
binding  of  LF  to  the  PA63  heptameric  bundle  by  incor¬ 
porating  multiple  copies  of  PI  and  P2  bound  to  a 
molecular  spacer  to  form  a  polyvalent  inhibitor.  These 
polyvalent  inhibitors  have  been  shown  to  provide  both 
in  vivo  (Fisher  344  rats,  75  nM  complete  protection)  and 
in  vitro  (IC50  approximately  20  nM,  CHO  cells)  protec¬ 
tion  against  anthrax  intoxication.2  In  our  laboratory  we 
are  interested  in  determining  the  conformation  adopted 
by  PI  on  binding  to  the  PA63  heptamer  bundle  in 
aqueous  solution  in  order  to  understand  the  structural 
and  physicochemical  requirements  for  binding.  This 
information  will  be  used  to  design  small  molecule 


Table  1.  Partial  Chemical  Shift  Assignments  (ppm)  for  PAS1 
(1.75  mg)  Bound  to  3.0  mg  of  PA63  (Molar  Ratio  25:1  Peptide  to 
Protein)  pH  =  4.0,  150  mM  Sodium  Acetate  Buffer, 
Temperature  =  277  K 


residue 

NH  oH 

PH 

yh 

other 

1  His 

3.942 

2.911 

Ar:  0.089/6.879 

2  Thr 

8.454  4.004  3.798 

0.792 

3  Ser 

8.297  4.111 

3.461/3.404 

4  Thr 

7.860  3.709  3.568 

0.591 

5  Tyr 

7.645  3.585  2.048/1.734 

Ar:  6.433/6.322 

6  Trp 

7.201  3.951 

2.899/2.803 

Ar  NH:  10.083" 

Ar:  6.868/6.796/6.962'’ 

7  Trp 

6.430  3.942 

2.598/1.909 

Ar  NH:  9.905" 

Ar:  6.614/6.768/6.458'’ 

8  Leu 

6.967  3.830 

1.004 

0.624 

0.390 

9  Asp 

7.651  4.125 

2.402 

2.284 

10  Gly 

7.890  3.486 

11  Ala 

7.720  4.078 

0.005 

12  Pro 

3.879 

1.818/1.555 

1.594/1.478 

3.254/3.167 

a  These  assignments  may  be  switched.  b  These  assignments  may 
be  switched. 


inhibitors  with  improved  therapeutic  index  compared 
to  PI.  Collier  and  co-workers  reported  that  at  basic  and 
neutral  pH,  PA63  forms  water  soluble  heptamer  bundles, 
and  at  acidic  pH  of  5.0  or  less,  PA6.3  forms  water  soluble 
heptamer  bundles  that  spontaneously  insert  into  lipid 
membranes.  It  was  therefore  assumed  that,  under  the 
conditions  employed  in  this  study,  the  majority  of  the 
PA63  is  involved  in  the  heptamer  bundle.11’12  No  ad¬ 
ditional  experiments  were  conducted  to  confirm  this 
assumption. 

Results 

The  first  step  in  this  investigation  was  to  determine 
whether  the  binding  of  PI  with  PA^a  occurs  under  the 
conditions  of  fast  or  slow  exchange.  The  amide  and 
aromatic  regions  of  the  !H  spectra  of  (A)  4.76  x  10  ~5  M 
PAe3,  (B)  4.76  x  10"5  M  PAea  plus  5.23  x  10  '4  M  PI 


PAe.r  Bound  Conformation  of  PI 
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Figure  2.  Top:  The  amide  to  amide  region  of  100  ms  mixing  time  NOESY  of  PI  in  buffer  at  277  K.  Bottom:  The  amide  to  amide 
region  of  the  100  ms  mixing  time  trNOESY  spectrum  of  1.75  mg  of  PI  (1.23  x  10-3  M)  with  3.0  mg  of  PA63  (4.76  x  10  5  M) 
(peptide  to  protein  ration  25:1)  in  1.0  mL  (90/10  H2O/D2O  150  mM  sodium  acetate  buffer,  pH  =  4.0).  The  one  letter  amino  acid 
symbols  are  used  to  indicate  amide  proton  to  amide  proton  interactions.  The  spectrum  was  collected  at  a  temperature  of  277  K. 
Note  that  key  i  to  i  4-  1  and  i  to  i  4-  2  NH-NH  connectivities  are  not  observed  in  the  trNOESY  spectrum  of  PI  in  buffer.  Thus 
the  NOEs  observed  in  the  trNOESY  spectrum  of  PI  bound  to  PA63  must  be  derived  from  the  transferred  NOE  effect. 


(peptide  to  protein  ratio  11:1),  (C)  4.76  x  10~5  M  PA63 
plus  1.23  x  10~3  M  PI,  (peptide  to  protein  ratio  25:1), 
and  (D)  8.3  x  10-4  M  PI  (all  in  150  mM  sodium  acetate 
buffer,  pH  =  4.0,  277  K)  are  given  in  Figure  1.  As  can 
be  seen  in  this  figure,  there  are  only  minor  detectable 
changes  in  chemical  shift  with  increasing  peptide 


concentration.  A  comparison  of  spectrum  C  with  the 
spectrum  of  PI  in  buffer  (spectrum  D)  shows  broadening 
of  the  resonance.  This  is  due  to  the  different  7*2 
relaxation  rates  of  the  peptide  bound  to  the  PA63 
compared  to  the  free  peptide.13  There  again  is  no 
detectable  change  in  amide  chemical  shifts.  This  obser- 


Figure  3.  Top:  The  amide  to  alkyl  region  of  50  ms  mixing 
time  NOESY  of  PI  in  buffer  at  277  K.  Bottom:  The  amide  to 
amide  region  of  the  50  ms  mixing  time  trNOESY  spectrum  of 
1.75  mg  of  PI  (1.23  x  10-3  mM)  with  3.0  mg  of  PA63  (4.76  x 
10“5  M)  (peptide  to  protein  ratio  25:1)  in  1.0  mL  (90/10  H20/ 
D20  150  mM  sodium  acetate  buffer,  pH  =  4.0).  The  one  letter 
amino  acid  symbols  are  used  to  indicate  amide  proton  to  amide 
proton  interactions.  The  spectrum  was  collected  at  a  temper¬ 
ature  of  277  K.  Note  that  key  i  to  i  +  1,  i  to  i  4-  2  and  i  to  i  + 
3  NH-aH  connectivities  are  not  observed  in  the  trNOESY 
spectrum  of  PI  in  buffer.  Thus  the  NOEs  observed  in  the 
trNOESY  spectrum  of  PI  bound  to  PAe3  must  be  derived  from 
the  transferred  NOE  effect. 


vation  is  consistent  with  “fast  exchange  behavior”.14 
Thus,  the  observed  NMR  signals  will  be  the  weighted 
average  of  the  contributions  of  the  PA63-bound  and  free 
conformers.15  Since  the  concentration  of  the  “free” 
peptide  is  approximately  one  and  one-half  orders  of 
magnitude  higher  than  the  protein-bound  peptide,  the 
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chemical  shift  values  observed  in  the  NMR  spectrum 
are  solely  representative  of  the  “free”  conformer.15 
However,  the  transferred  nuclear  Overhauser  effect  is 
very  sensitive  to  NOEs  induced  onto  the  free  ligand  by 
binding  to  the  larger  protein.16-18  These  NOEs  are 
observed  in  the  presence  of  “fast  exchange”,  and  they 
may  be  employed  to  develop  information  concerning  the 
protein-bound  conformation  of  PI.1419 

The  second  step  in  this  investigation  was  to  determine 
the  protein-bound  conformation  of  PI.  Using  the  2D- 
TOSCY20-22  and  2D-NOESY23  spectra,  the  !H  chemical 
shifts  for  PI  (Table  1)  were  obtained  by  employing  the 
sequence-specific  resonance  assignment  technique  de¬ 
veloped  by  Wiithrich  and  co-workers. 24  The  PA63-bound 
conformation  of  PI  was  determined  using  five  trNOESY 
experiments.  The  magnetization  of  smaller  molecules, 
like  PI  (MW  1433),  on  binding  to  large  biomolecules 
such  as  PA63  (MW  63000)  is  labeled  by  the  magnetiza¬ 
tion  of  the  large  biomolecules.  As  a  result,  these  smaller 
molecules  yield  large,  negative  NOEs  that  build  up 
rapidly  and  reach  a  maximum  intensity  at  relatively 
short  mixing  times.  Thus,  the  binding  of  a  small 
molecule  to  a  large  molecule  can  be  detected  by  the  “size, 
sign  and  build-up  rates”14  of  NOEs  in  the  trNOESY 
spectrum.  Negative  peaks  were  observed  in  the  tr¬ 
NOESY  spectra  of  PI  in  buffer  and  in  the  presence  of 
PA(53,  but  in  the  latter  case  the  peaks  were  more  intense. 
Furthermore,  as  illustrated  in  Figures  2,  3,  and  4,  the 
trNOESY  spectrum  of  PI  bound  to  PA63  shows  many 
more  NOE  cross-peaks  than  are  observed  in  the  tr¬ 
NOESY  spectrum  of  PI  in  buffer  under  the  identical 
conditions.  In  addition,  the  NOE  build-up  rates  for  PI 
bound  to  PA63  are  very  different  from  the  NOE  build¬ 
up  rates  of  PI  in  buffer.  The  NOE  build-up  rates  for 
seven  intraresidue  and  interresidue  proton  pairs  of  PI 
bound  to  PAe3  and  three  of  these  proton  pairs  for  PI  in 
buffer  are  given  in  Figure  4.  For  PI  bound  to  PA63,  all 
NOE  intensities  reach  a  maximum  value  in  the  range 
of  50  to  100  ms,  while  those  for  PI  in  buffer  reach  their 
maximum  intensities  in  the  range  of  300  to  400  ms. 
These  observations  indicate  that  the  NOEs  observed  for 
PI  bound  to  PA(33  are  transferred  NOEs. 

The  following  NOE  connectivities  in  the  trNOESY 
spectra  of  this  PI  bound  to  PA63  are  sufficient  to 
characterize  a  helical  secondary  structure  involving 
residues  3-9  of  the  peptide:  (1)  four  sequential  medium 
NH-NH/+1  (residues  2—6)  and  four  sequential  medium 
NH-NH/+1  (residues  7-11)  (2)  a  series  of  sequential 
medium  NH-NH/+2  (residues  2-4,  5-7,  6-8,  7-9, 
8-10),  (3)  eleven  sequential  CaH-NHi'+l,  (4)  four  CaH- 
NH/+2  (residues  1-3,  4-6,  5-7,  9-11),  and  (5)  three 
sequential  medium  CaH-NHi+3  (residues  4-7,  5-8, 
6-9).  These  NOEs  are  given  in  Figure  6. 

The  NOEs  observed  in  the  trNOESY  spectra  of  PI 
bound  to  PA(53  were  used  as  restraints  in  NOE  re¬ 
strained  simulated  annealing  calculations  using  a  pre¬ 
viously  reported  protocol  to  generate  100  structures  of 
PI  bound  to  PA63.25  Thirty  of  the  lowest  energy  struc¬ 
tures  (average  energy  350.8031  kcal/mol)  are  shown  in 
Figure  7  superimposing  the  backbone  atoms  of  residues 
3-10  onto  each  other.  The  rmsd  for  the  superimposi¬ 
tions  of  the  backbone  atoms  of  residues  3-10  ranges 
from  0.206  to  1.35  A  (0.841  A  average).  The  average 
dihedral  angles  for  these  structures  are  given  in  Table 
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Y5-W7 


L8-W7 


W6-W7 


'  '  '  ~ '  '  n - - - 1 — * . i  - -i . i 

7-5  7-°  6-5  PPn  7.5  7.0  6.5  ppm 

Figure  4.  Scaled  slices  taken  from  the  100  ms  trNOESY  corresponding  to  the  chemical  shift  of  the  Trp7  NH  resonance  (left)  PI 
bound  to  PA«:3  and  (right)  PI  in  buffer  only.  This  data  clearly  indicates  that  strong  NH-NH  NOEs  are  observed  in  the  presence 
of  P Ae.-j  only. 


Table  2.  Average  Dihedral  Angles  for  the  30  Low  Energy 
Conformations  for  PI  Bound  to  the  PA6.3  Heptamer  Bundle 


residue 

av  (j> 

av  ip 

residue 

av  <p 

av  xp 

Hisl 

57.28 

Trp7 

-67.29 

-40.38 

Thr2 

-72.71 

-27.03 

Leu8 

-73.54 

-48.07 

Ser3 

-75.18 

-29.33 

Asp9 

-116.07 

29.60 

Thr4 

-79.80 

-20.19 

GlylO 

-14.11 

-11.65 

Tyr5 

-75.18 

-39.89 

Alall 

-89.30 

130.99 

Trp6 

-70.35 

-34.45 

Pro  12 

2.  Clearly,  the  midregion  of  the  peptide  consisting  of 
residues  4-8  converges  to  an  a-helix.  Residues  3  and 
9,  while  maintaining  a  helical  conformation,  appear  to 
be  somewhat  more  flexible.  The  remaining  N-terminal 
and  C-terminal  residues  appear  to  be  very  flexible. 

Discussion 

Our  NMR  and  NOE  restrained  simulated  annealing 
calculations  indicate  that  PI  adopts  a  helical  structure 
involving  residues  3-9  on  binding  to  the  PA63  heptamer 
bundle.  Collier  and  co-workers  have  suggested  that  the 
tetrapeptide  sequence  Try5-Trp6-Trp7-Leu8  plays  an 
important  role  in  binding  to  a  hydrophobic  pocket  near 
the  LF  binding  site  of  PA63.  Our  results  indicate  that 
the  tetrapeptide  sequence  exists  in  the  bound  conforma¬ 
tion  as  a  very  well  defined  a-helical  structure  (Figure 
7).  These  four  residues  are  on  different  faces  of  the  helix 
and  are  not  interacting  with  PAe3  in  a  linear  fashion. 
Therefore,  it  must  be  concluded  that  the  relative  posi¬ 
tions  of  the  side  chains  in  three-dimensional  space  are 
of  critical  importance  in  binding  to  PAe3.  On  the  basis 
of  examination  of  this  structure,  we  feel  that  the 
N-terminal  His1  residue  (which  is  also  common  to  PI 


and  P2),  even  though  it  is  more  conformationally  flexible 
than  residues  3—9,  plays  a  critical  role  in  binding  by 
providing  an  electrostatic  interaction  on  the  same  face 
of  the  helix  as  Tyr5  and  Trp6  residues.  In  fact,  as  can 
be  seen  in  Figure  8,  the  side  chains  of  all  three  residues 
are  very  close  to  each  other  (within  7  A).  This  is  in 
contradiction  to  the  structure  proposed  by  Glick  and  co¬ 
workers.26  This  group  investigated  the  binding  interac¬ 
tion  of  the  tetrapeptide  Try5-Trp6-Trp7-Leu8  using  con¬ 
jugate  gradient  minimization  computational  methods 
with  the  crystal  structure  of  the  PA63  heptamer  “pre- 
pore”  reported  by  Petosa  and  co-workers.11  Their  results 
indicate  that  the  tetrapeptide  binds  in  an  extended 
^-conformation  to  a  hydrophobic  pocket  on  PA^-formed 
residues  Trp226,  Tyr462,  and  Phe464.  This  discrepancy  is 
not  unexpected,  since  the  computational  methods  did 
not  involve  the  complete  amino  acid  sequence  of  the 
peptide,  particularly  the  N-terminal  His1  and  the  Asp9 
residues.  The  electrostatic  potential  map  of  the  PA43 
heptamer  bundle  bound  conformation  of  PI  is  shown 
in  Figure  9.  Blue  regions  represent  positive  electrostatic 
potential  associated  with  the  His1  residue;  the  red 
regions  represent  negative  electrostatic  potential  as¬ 
sociated  with  the  Asp9  residue.  The  tetrapeptide  se¬ 
quence  Tyr5-Trp6-Trp7-Leu8  connecting  these  two  re¬ 
gions  is  electrically  neutral.  Since  the  exact  binding  site 
of  PI  is  still  unknown,  it  may  be  reasonable  to  speculate 
from  this  electrostatic  potential  profile  that  the  comple¬ 
mentary  binding  site  requires  a  large  distribution  of 
hydrophobic  regions  (white  portion,  Figure  9)  flanked 
by  a  negative  potential  region  (red  regions)  for  favorable 
interaction  with  the  His1  residue  and  a  positive  poten- 
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Figure  5.  Top:  NOE  build-up  rates  for  seven  intraresidue 
and  interresidue  proton  pairs  of  PI  bound  to  PA6.3.  (Dark  blue 
triangles:  Thr}NH-Ser\iH.  Green  triangles:  Se^NH-Thr^aH. 
Light  blue  tilted  squares:  Thr4NH-aH.  Black  squares:  Thr2- 
NH— aH.  Red  circles:  Ser3NH-aH.  Magenta  triangles:  Leu8- 
NH-aH.  Yellow  triangles:  Leu8N-Tyr5aH.)  The  NOEs  reach 
a  maximum  negative  intensity  in  the  mixing  time  range  of  50 
to  100  ms.  This  behavior  is  consistent  with  that  observed  with 
the  transferred  NOE  effect.14  Bottom:  NOE  build-up  rates  for 
three  intraresidue  and  interresidue  proton  pairs  of  PI  in  buffer 
only.  (Red  circles:  Ser3NH— Thr2aH.  Black  squares:  Thr2NH— 
His^H.  Green  triangles:  Ser3NH— aH.)  The  NOEs  reach  a 
maximum  negative  intensity  in  a  mixing  time  range  of  300  to 
400  ms.  This  behavior  is  consistent  with  a  “normal”  NOE 
effect.  Note  that  NOE  intensities  for  the  two  graphs  are  not 
scaled  relative  to  each  other. 

tial  region  (blue)  for  the  interaction  with  Asp9  residue. 
Thus,  the  electrostatic  potential  map  clearly  demon¬ 
strates  that  any  attempt  to  dock  only  the  tetrapeptide 
to  PA(53  will  not  represent  the  complete  physicochemical 
surface  of  PI. 

The  overall  structural  requirements  for  binding  of  PI 
with  the  PA63  heptamer  bundle  have  yet  to  be  investi¬ 
gated.  The  importance  of  the  common  hydrophobic 
sequence  YWWL  for  binding  to  PA6.3  will  be  determined 
using  selective  amino  acid  replacement.  In  our  opinion, 
the  hydrophobic  character  of  the  side  chains  of  the  Tyr 
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Figure  6.  Key  NOE  connectivities  used  as  interproton 
distance  restraints  in  the  simulated  annealing  calculations  to 
determine  the  conformation  of  peptide  1  bound  to  the  PA«3. 
Line  thickness  indicates  the  strength  of  the  NOE.  Strong, 
medium,  and  weak  correspond  to  an  upper  distance  of  strong, 
2.8  A,  medium  3.5  A,  and  weak  to  very  weak  4.0  and  5.0  A.34 

and  Trp  residues  and  the  ability  to  act  as  H-bond 
acceptor/donator  must  be  systematically  investigated. 
Therefore  each  Tyr  and  Trp  residue  will  be  replaced  in 
turn  with  Phe.  Because  the  hydrophobic  volume  of  Trp 
is  larger  than  Phe,  each  Trp  residue  will  also  be  replaced 
with  a  naphthylalanine  residue  (Nal).  Our  hypothesis 
concerning  the  importance  of  potential  electrostatic 
interactions  of  the  N-terminal  His  residue  with  the  PA63 
heptamer  bundle  will  be  tested  by  synthesis  of  two 
analogues  of  PI.  The  His  residue  will  be  replaced  with 
the  neutral  residue  Gly,  and  by  the  acidic  residue  Asp. 
Therefore  to  better  define  the  structure  function  of  PI 
we  proposed  to  synthesize  the  analogues  given  in  Table 
3.  NMR  and  molecular  modeling  techniques  will  be  used 
to  determine  the  effect  of  modification  on  the  PA&3- 
bound  conformation  of  each  analogue. 

Materials  and  Methods 

Materials.  PA63  was  purchased  from  List  Biological  Labo¬ 
ratories  (Campbell,  CA).  PI  was  prepared  using  Fmoc-based 
solid  phase  peptide  synthesis  by  Biopeptide  Company  (San 
Diego,  CA).  Deuterated  sodium  acetate  and  deuterated  acetic 
acid  were  purchased  from  Sigma  Chemical  Company  (St. 
Louis,  MO).  Isotopically  enriched  2H20  was  purchased  from 
Aldrich  Chemical  Co.,  Inc.  (St.  Louis,  MO).  All  materials  were 
used  without  further  purification. 

Sample  Preparation.  Samples  were  prepared  for  the  NMR 
binding  studies  using  1.0  mL  aliquots  of  solutions  composed 
of  3.0  mg  (4.76  x  10  5  M)  of  PAe3  and  1.7  mg  (1.23  x  10  3  M) 
of  PI,  in  9:1  1H20/2H20  buffered  with  150  mM  sodium  acetate 
to  a  pH  of  4.0.  For  the  reference  studies  a  solution  of  1.7  mg 
(1.23  x  10  3  M)  of  PI  in  9:1  ^O^HiO  buffered  with  150  mM 
sodium  acetate  to  a  pH  of  4.0  was  prepared. 

NMR  Experiments.  All  *H  NMR  data  was  collected  using 
a  Bruker  Avance-600  NMR  spectrometer.  The  application  of 
2D-NMR  methods  to  investigate  interactions  between  enzymes 
and  substrates  via  the  trNOESY  experiment  has  been  exten¬ 
sively  reviewed  in  the  literature  and  will  not  be  dis¬ 
cussed.131416- 1927-31  The  proton  spectrum  of  PI  was  assigned 
by  the  interactive  analysis  of  the  2D  TOCSY2022  and  tr¬ 
NOESY14’23  spectra.  The  TOCSY  experiment  employed  a 
modified  MLEV-1721  spin-lock  sequence  for  a  total  mixing  time 
of  80  ms,  including  the  2.5  ms  trim  pulses  at  the  beginning 
and  the  end  of  the  spin-lock  sequence.  Five  trNOESY  spectra 
of  PI  bound  to  PAe.3  were  collected  with  mixing  times  of  10, 
25,  50,  100,  and  200  ms.  Five  NOESY  spectra  of  PI  in  buffer 
were  collected  with  mixing  times  of  50,  100,  200,  300,  and  400 
ms.  Solvent  suppression  in  all  experiments  was  accomplished 
by  application  of  the  WATERGATE  (WATER  suppression  by 
GrAdient  Tailored  Excitation)  pulse  sequence  developed  by 


PAes-Bound  Conformation  of  PI 


Journal  of  Medicinal  Chemistry ,  2004 ,  Vol.  47,  No.  22  5353 


Table  3.  Analogues  of  PI  To  Be  Prepared  To  Evaluate  the  Contributions  of  the  YWWL  Tetrapeptide  Sequence  and  the  His  Residue 
to  Binding  Affinity  to  the  PA63  Heptamer  Bundle 


amino  acid  sequence 


function 


His-Thr-Ser-Thr-Tyr-Trp-Trp-Leu-Asp-Gly-Ala-Pro 

His-Thr-Ser-Thr-Phe-Trp-Trp-Leu-Asp-Gly-Ala-Pro 

His-Thr-Ser-Thr-Tyr-Phe-Trp-Leu-Asp-Gly-Ala-Pro 

His-Thr-Ser-Thr-Tyr-Trp-Phe-Leu-Asp-Gly-Ala-Pro 

His-Thr-Ser-Thr-Tyr-Trp-Trp-Phe-Asp-Gly-Ala-Pro 

His-Thr-Ser-Thr-Tyr-Nal-Trp-Leu-Asp-Gly-Ala-Pro 

His-Thr-Ser-Thr-Tyr-Trp-Nal-Leu-Asp-Gly-Ala-Pro 

Ala-Thr-Ser-Thr-Tyr-Trp-Trp-Leu-Asp-Gly-Ala-Pro 

Asp-Thr-Ser-Thr-Tyr-Trp-Trp-Leu-Asp-Gly-Ala-Pro 


PI 

H-bonding 

H-bonding 

H-bonding 

hydrophobic  effect  ( Phe  vs  Leu ) 
H-bonding  and  bulk 
H-bonding  and  bulk 
electrostatic  interaction:  neutral 
electrostatic  interaction:  negative 


Figure  7.  Stereoview  of  the  30  lowest  energy  structures  (average  energy  350.8031  kcal/mol)  of  the  PA63-bound  conformation  of 
PI  superimposing  the  backbone  atoms  of  residues  3  to  10  onto  each  other.  The  rmsd  for  this  superimposition  ranges  from  0.206 
to  1.35  A  (0.841  A,  average).  Clearly  the  midregion  of  the  peptide  consisting  of  residues  4-8  converges  to  an  a-helix.  Residues  3 
and  9  of  the  helix  appear  to  be  somewhat  more  flexible.  The  remaining  N-terminal  and  C-terminal  residues  appear  to  be  very 
flexible. 


Figure  8.  A  stereoview  of  a  representative  structure  of  PI  bound  to  the  PA^s  heptamer  bundle.  Only  heavy  atoms  are  shown. 
The  tetrapeptide  Try5-Trp6-Trp7-Leu8,  believed  to  bind  with  PAe.3  via  hydrophobic  interactions,  is  shown  in  dark  blue.  The  N-terminal 
His1  residue,  which  we  feel  plays  a  critical  role  in  electrostatic  interactions,  is  shown  in  red.  Intramolecular  distances  between 
residues  are  given  to  highlight  the  closeness  of  particular  residues  to  each  other. 


Sklenar  and  co-workers.32  2D-TOCSY  spectra  were  collected 
with  64  transients  and  1024  t\  increments  with  a  1.5  s  recycle 
delay  for  all  experiments.  2D-trNOESY  spectra  were  collected 
with  128  transients  and  1024  t\  increments.  Data  were 
collected  at  277  K.  The  spectral  width  was  9090.9  Hz  in  both 


dimensions.  All  spectra  were  acquired  with  2048  data  points 
in  F2.  All  chemical  shifts  were  referenced  internally  to  DSS 
(0.00  ppm).  Spectra  were  processed  using  XWINNMR  ( Bruker) 
on  a  Silicon  Graphics  O2  workstation.  Each  data  set  was 
multiplied  by  a  Gaussain  window  function  with  —10  Hz  line 


5354  Journal  of  Medicinal  Chemistry ,  2004,  Vol.  47,  No.  22 


Hicks  et  al. 


Figure  9.  Electrostatic  potential  map  of  the  PA63  heptamer 
bundle  bound  conformation  of  PI.  Blue  regions  represent 
positive  electrostatic  potential  associated  with  the  His1  residue; 
the  red  regions  represent  negative  electrostatic  potential 
associated  with  the  Asp9  residue.  The  tetrapeptide  sequence 
Tyr5-Trp6-Trp7-Leu8  connecting  these  two  regions  is  electro¬ 
statically  neutral.  This  map  clearly  shows  that  any  attempt 
to  dock  only  the  tetrapeptide  to  PA63  will  not  represent  the 
physicochemical  surface  of  PI. 


broadening  in  each  dimension  before  transformation  to  pro¬ 
duce  matrices  consisting  of  2048  data  points  in  both  dimen¬ 
sions. 

Molecular  Modeling.  All  molecular  modeling  experiments 
were  performed  using  the  Biopolymer,  NMR  Refine  Advanced, 
and  Discover  modules  of  the  Insightll  package  from  Accelrys,33 
operating  on  a  Silicon  Graphics  Octane  workstation.  Data  from 
the  five  WATERGATE-trNOESY  experiments  were  used  to 
classify  peak  volumes  over  a  range  of  strong  to  very  weak, 
corresponding  to  the  upper-bound  interproton  distance  re¬ 
straints  of  2.8,  3.5,  4.0,  and  5.0  A  29,34-36  with  strong  NOEs 
(2.8  A)  taken  from  the  10  and  25  ms  WATERGATE-trNOESY, 
medium  (3.5  A)  from  the  50  and  100  ms  WATERGATE- 
trNOESY,  and  weak  (4.0  and  5.0  A)  from  the  200  ms  WA¬ 
TERGATE-trNOESY,  respectively.  Since  stereospecific  assign¬ 
ments  could  not  be  made  initially,  pseudoatoms  were  employed 
using  the  center-of-mass  approach,  with  both  intraresidue  and 
long-range  correction  factors  added  to  the  distance  restraints.25 
In  the  case  of  terminal  methyl  groups,  a  distance  correction 
of  0.5  A  was  added  to  the  upper  limits  for  distance  restraints.35 
After  one  set  of  MD/SA  calculations,  stereospecific  assignments 
were  made  by  manually  measuring  interproton  distances  and 
selecting  those  which  were  within  the  acceptable  distance 
range.  Correction  factors  were  removed,  and  MD/SA  calcula¬ 
tions  were  repeated.  A  final  set  of  100  structures  were 
generated  for  the  PA63  heptamer  bound  PI  starting  from 
templates  with  completely  randomized  coordinates  followed 
by  the  application  of  a  MD/SA  protocol  similar  to  that  used 
by  us  earlier,  for  determining  the  SDS  bound  conformations 
of  members  of  the  tachykinin  family  of  neuropeptides.25,35 

Molecular  Electrostatic  Potentials  (MEP).  Molecular 
electrostatic  potential  (MEP)  calculations  on  the  peptides  were 
performed  using  the  DelPhi  module  of  Insightll.33  The  algo¬ 
rithm  allows  calculation  of  electrostatic  potential  in  and 
around  a  protein  to  be  displayed  as  contours  to  gain  qualitative 
information  on  protein -substrate  interactions.37-40  The  gener¬ 
ated  MEP  profile  may  be  considered  as  an  “interaction 
pharmacophore”  of  the  molecule.40 
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